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Abstract. We present a novel theoretical tool to analyze the dynamical behaviour of a Be disk 
fed by non-constant decretion rates. It is mainly based on the computer code HDUST, a fully 
three-dimensional radiative transfer code that has been successfully applied to study several 
Be systems so far, and the SINGLEBE code that solves the ID viscous diffusion problem. We 
have computed models of the temporal evolution of different types of Be star disks for different 
dynamical scenarios. By showing the behaviour of a large number of observables (interferometry, 
polarization, photometry and spectral line profiles), we show how it is possible to infer from 
observations some key dynamical parameters of the disk. 
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1. Presentation of the simulations 

In order to analyze the observational signatures of different dynamical scenarios, we use 
the SINGLEBE code ( [Okazaki, 2007[ ) that computes the temporal evolution of the surface 



density for a given stellar mass loss history and an a viscosity parameter ( Shakura & Sunyaev, 1973 ) 
of the disk. We then use the surface density at a given time as an input for a three- 
dimensional non- LTE Monte Carlo code called HDUST QCarciofi k Bjorkman, 20061 ). 
By repeating this for different epochs of the disk evolution, we can follow the evolution 
of several observables (SED, polarization, images, etc). 

A wide range of observables can be derived from the HDUST simulations. By means 
of a comparative work, they can unveil the observational features of a specific dynamical 
scenario occurring in a viscous disk. The final objective of these quantities is to be an- 
alyzed and compared to real polarimetric, photometric, spectroscopic or interferometric 
observations. This method has the potential to allow one observer to unambiguously infer 
the size, the a parameter and the mass loss history of the observed system. 



2. Some results 

As preliminary results, we present some observables and their correlation which clearly 
show different temporal behaviour depending on the a parameter (FiglTJ) . We present 
results for a dynamical scenario in which 3 year long outbursts are separated by 3 year 
long quiescent phases. The influence of a is understood by considering that the larger a, 
the faster the viscous diffusion occurs. 
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Figure 1. Upper panel: polarization change across the Balmer Jump (BJ) vs V band polarization 
for two a parameter values (0.1 in dashed line, 0.7 in full line) and at an inclination angle of 70 
degrees. Epochs are marked and counted in years, the correlation loops move clockwise. Lower 
panel: temporal evolution of these two quantities. The thickest line shows the mass loss history 
arbitrarily scaled to the range of the graphic. It has been turned on and off every 3 years from 
to 39 years. That figure depicts the 32-39 years period. 

3. Conclusion and Perspectives 

These simulations can predict the observational signatures of dynamical scenarios by 
generating a wide range of observables. Some of these observables are thus powerful tools 
to estimate some key parameters of observed Be systems such as the a viscosity param- 
eter. This work is currently ongoing and more complex dynamical scenarios have to be 
investigated ( |Haubois et al., in preparation! ) . Moreover, the results of these studies being 
more easily understandable with videos, we plan to make them accessible to the commu- 
nity in a near future through a dedicated website ([http://www.astro.iag.usp.br/beacon7). 
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